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ABSTRACT 

The main objective of this paper is to study the possibility of triggered star formation on the 
border of the H ll region S233, which is formed by a B-star. Using high-resolution spectra we 
determine the spectral class of the ionizing star as BO.5 V and the radial velocity of the star to 
be —17.5 ±1.4 km s“^. This value is consistent with the velocity of gas in a wide field across 
the S233 region, suggesting that the ionizing star was formed from a parent cloud belonging 
to the S233 region. By studying spatial-kinematic structure of the molecular cloud in the S233 
region, we detected an isolated clump of gas producing CO emission red-shifted relative to 
the parent cloud. In the UKIDSS and WISE images, the clump of gas coincides with the 
infrared source containing a compact object and bright-rimmed structure. The bright-rimmed 
structure is perpendicular to the direction of the ionizing star. The compact source coincides in 
position with IRAS source 05351-b3549. All these features indicate a possibility of triggering 
formation of a next-generation star in the S233 region. Within the framework of a theoretical 
one-dimensional model we conclude that the ‘collect-and-collapse’ process is not likely to 
take place in the S233 region. The presence of the bright-rimmed structure and the compact 
infrared source suggest that the ‘collapse of the pre-existing clump’ process is taking place. 
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1 INTRODUCTION 

A considerable number of stars in our and other galaxies are formed 
as a result of star formation triggered by the expansion of H II 
regions. Such expanding H II regions have a great variety. They 
can appear as a result of rare energeti c events like a sup ernova ex¬ 
plosion producing huge superbubbles (iKang et al.ll20121) , rare ener¬ 
getic objects like WR stars o r clusters of massive stars (lOev et al.l 
I2OO5I : I Schneider et al.|[201^ . These energetic objects attract the 
majority of effort in the field of triggered star formation. At the 
same time star formation can be efficient on the bo rders of the H II 
regions formed by almost isolated massive s tars (iKirsanova et al.l 
I2OO8I : ISnider et ^l2009l : IZavagno et al.ll201Qh . 

Rather frequently triggered star formation is realized as a 
collapse of a pre-existing clump. In this process contraction of 
a pre-existing molecular clumps is initiated by the photoioniza¬ 
tion and shock fronts propagating away from the H II region. 
Bright rims are expected to appear at the sides of the clumps 
which face the ionizing star (simila r to those observed e.g. by 
iThompson. Urquhart. & Whit3 (l2004lfi . Clear examples of the star 
formation due to compression of pre-existing dense condensa- 
tions are found in IC 1396, IC 1805, SFO 79 star-forining re gions 

jWeikard et al.lll99d : iHever et ^1 19961 : lUrauhart et al.ll2004h . 

Another type of triggered star formation is the ‘collect-and- 


collapse’ process, first proposed bv lElmegreen & Lad3 (Il977h . In 
this process massive OB-stars create an H II region expanding into 
the ambient molecular cloud due to the presence of a pressure dif¬ 
ference between the molecular and ionized gas. This expansion 
compresses ambient molecular material, consecutive fragmentation 
of which creates molecular clumps and filaments. These clumps are 
progenitors for the next generations of stars. The important aspect 
of this scenario is that there is a delay between the expa nsion of the 
H II region and the onset of triggered star formation jEImegreenl 
Il998h . This is different from the case of star formation by collapse 
of pre-existing clumps where the formation of the star starts almost 
immediately. Clear examples of star formation by the ‘collect-and- 
collapse’ process in our Gala xy can be found in SI04 and S212 
(jPeharveng et al.l200ll2008lj . 

I Whit worth et al.l(l 19941) investigated properties of the ‘collect- 
and-collapse’ process of star formation with a one-dimensional 
model of the shocked shell driven by expansion of the H II region 
into a uniform medium. They concluded that resulting fragments 
are likely to have high mass (> 7 M©). It was found that gravita¬ 
tional fragmentation of the shell takes place when the c olumn den¬ 
sity of hydrogen in the shell reaches 6 X 10^^ cm ^. IPale et^ 
(l2009h use the thee-dimensional smooth particle hydrodynamic 
(SPH) model to test the applicability of the thin-shell approxima- 
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tion used in i Whitworth et al.l(ll994) for more realistic astrophysical 
situation. They concluded that this approximation should be used 
with some caution because of the problem with the boundary con¬ 
ditions which greatly influence the fragmentation of the shell. 

In this paper, we consider the border of an H II region around 
an early-type star. It will be shown that there are signs of the con¬ 
siderable influence of this ionizing star on the ambient molecular 
cloud. This influence possibly results in triggering the formation of 
the next-generation star(s) on the border of the S233 H II region. 

The S233 region is a small 2' in diameter) optical 
emission nebula formed by the ionizing star USNO-A2 1200- 
03588518 (aj 2 ooo = 05'’38"’31!5, 5j2ooo = +35°51'19") with 
visual magnitude V = 11.7"^. The S233 region is a part of 
S231-235 star formation complex. This complex is located in the 
Perseus spiral arm and contains four devel oped H II regions: S23I, 
S232, S233, and S235 (iHever et al.lll996h . Published distance es- 
timates to the objects in this complex vary from 1.6 to 2.3 kpc 
teeipurth & YanH2008 ). The S2 33 region is mentioned in the re¬ 


view of iReipurth & YanI (l2008h . which states that the S2 33 H II 
region is excited by a B 1.5 II star (iHunter & Masseviri99Clh . In the 
present study we conducted a spectroscopic analysis of the ionizing 
star using high-resolution spectroscopic data and have shown that 
it is a main sequence star of the B0.5 spectral type with mass of 


13 z b 1 Me. _ 

IWouterloot & Bran3 Il989h studied CO(l-O) molecular line 
emission toward the IRAS sources beyond the solar circle. They 
distinguished two separate components of the CO(l-O) line to¬ 
ward IRAS 05351-1-3549 situated within the optical S233 emis¬ 
sion. Our newly obtained data enable us to investigate the spatial 
structure of these CO components. The main gas component of 
CO emission (12.1 zb 0.4 K at —18.2 zb 0.4 km s“^) has a large- 
scale spatial extension, but the second component (2.6 zb 0.4 K at 
— 12.9 zb 0.3 km s“^) is a rather compact clump of gas situated 
close to the infrared source IRAS 05351-1-3549. 


2 OBSERVATIONS AND DATA REDUCTION 

2.1 The CO and CO Molecular Gas 

For our study of molecular gas in the vicinity of the S233 re¬ 
gion, we used information on emission in four spectral lines of the 
CO molecule: ^^CO(l-O), ^^CO(2-l) and the lines of isotopologue 
^^CO(l-O) and ^^CO(2-l). These lines are known to be an efficient 
probe of morphology and physical parameters of molecular gas. 

The maps of the ^^CO(2-l) and ^^CO(2-l) spectral lines were 
obtained with the 10-meter Heinrich Hertz Submillimeter Tele¬ 
scope (SMT), a facility of the Arizona Radio Observatory. The 
maps include a 70' x 50' region centered at RA(2000) = 5^40™, 
Dec.(2000) = 35° 50' which covers more than 20' around the 
vicinity of the S233 region. Observations were obtained in January 
2010. The whole map consists of subflelds with 10' x 10' size. 
Each subfleld was mapped using the On-The-Fly mode of the SMT 
telescope. The beam size was 32" for the ^^CO(2-l) and 33.5" 
for the ^^CO(2-l). Observed frequencies were 230.53800 GHz for 
the ^^CO(2-l) and 220.39868 GHz for the ^^CO(2-l) lines, re¬ 
spectively. The main beam efficiency was 0.81 for the ^^CO line 
and 0.78 for the ^^CO line. The observations were made with a 
10" angular spacing, which was at least three times smaller than the 
beam size. This allows avoiding sudden changes in intensity from 
pixel to pixel and positively affects the accuracy of determination 
of physical parameters in the molecular cloud. The spectral resolu¬ 
tion was 0.325 km for the ^^CO(2-l) line and 0.340 km 


and ^^CO(l-O) 


for the ^^CO(2-l) line. Achieved RMS noise levels were 0.23 K 
and 0.21 K for the ^^CO(2-l) and ^^CO(2-l) lines, respectively. 

The ^^CO(l-O) and ^^CO(l-O) spectral lines were mapped 
using the 13.7-meter FCRAO radio telescope with the 32-pixel 
SEQUOIA focal plane array receiver. The S231-S235 complex 
was mapped in January 2000. The map of each isotope line in¬ 
cludes a 150' X 150' region based on the galactic coordinate sys¬ 
tem, centered at I = 173.25°, b = 2.75° (RA = 5^40"^, 
Dec. = 36°07'). The main beam sizes were 45" and 47" for the 
^^CO(l-O) and ^^CO(l-O) lines, respectively. Observed frequencies 
were 115.27120 GHz and 110.20135 GHz for the ^^CO(l-O) and 
^^CO(l-O) lines, respectively. The pixel spacing was 22.5" and the 
velocity channel spacings were 0.127 km and 0.133 km 
for the ^^CO(l-O) and ^^CO(l-O) lines, respectively. Achieved 
RMS noise levels were 1.1 K and 0.63 K for the ^^CO(l-O) and 
^^CO(l-O) lines, respectively. 

The main beam efficiencies for the 
lines were 0.45 and 0.48, respectively. These efficiencies applied 
if the source of emission was approximately the same size as the 
main beam. However, if the emission is widely distributed on scales 
of 0.5° or more, it Alls the error beam of the FCRAO telescope, 
and in such a case one need to apply scattering efficiency of 0.7 
for both lines. Given the large range of angular sizes of the struc¬ 
tures in S233 complex neither main beam nor scattering efficiencies 
will give good calibration. To calibrate the FC RAO CO(l-O) data 
we used a deconvolution method described in (IPineda et al.ll2010h 
to remove a 0.5° error beam component and then we used main 
beam efficiencies. The deconvolution method involves a division 
of Fourier transform of source data by Fourier transform of the er¬ 
ror beam. Correction for the error beam improves the calibration of 
the line intensities by 10-30 percent. 

We used the _ MIRIAD software package 

jSault. Teuben. & WrightI Il995h for man ipulating data cubes 
and the KARMA package ( Gooch! Il995h for visualization. In 
this paper, we have extracted a 10' x 10' region centered at the 
position of the ionizing star in the S233 region with coordinates 
RA(2000) = 5^38""3R5, Dec.(2000) = 35°51'19". Because 
the data for different transitions was obtained with different 
receivers, we have applied some corrections to the original data for 
comparing line profiles observed with different telescopes. 

These corrections include the convolution of all data cubes 
to the same beam size of FWHM= 47". For the ^^CO (1-0) line 
this is the original beam size, thus no changes were made to this 
data cube, but for the ^^CO (2-1) and ^^CO (2-1) lines the original 
beam size is 32" and 33.5", respectively. Thus to achieve 47" beam 
size we convolved these line maps with ^47^ — 32^ = 34.4" and 
V 472 — 33.52 _ 32.9" Gaussians, respectively. For investigating 
the spatial-kinematic structure of molecular cloud the corrections 
are not necessary. 


2.2 Echelle-spectrum of the ionizing star 

The observations were obtained in January 2012 with the Echelle- 
spectrograph mounted on the Special Astro physical Observatory 
(SAO) 6-m telescope jPanchuk et al.l 1200^ . The detector of the 
Echelle-spectrograph was a large size 4608 x 2048-pixel CCD 
array with an image sheer. We used a setup yielding a resolving 
power of R = 50000, with a spectral coverage of 4000-7000 A. A 
hollow cathode Th-Ar lamp was used for the wavelength calibra¬ 
tion. 

We have obtained 6 spectra with a 45 minutes exposure time 
for each spectrum. Initial inspections showed no significant veloc- 
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Figure 1. Smoothed and rebinned to R=4000 spectrum of S233 ionizing 
star (lo wer spectrum) and spec t rum o f HD 36960 (BO.5 V) from the digital 
atlas of lWalborn & FitzpatrickI Il99(t) . 


ity variation in this sample. Thus, all spectra were combined to im¬ 
prove the signal-to-noise ratio (SNR) and a master spectrum was 
created for subsequent analyses. The initial steps of the data reduc¬ 
tion process (removal of cosmic ray features, background subtrac¬ 
tion, and spectral order extraction) we re done under the EC HELLE 
context of MIDAS software package jYushkin et aLll200^ . while 
the final steps (i.e., continuum normalization, radial velocity and 
equivalent widths (EW) measurements ) were performed using the 
package DECH20 (iGalazutdinovI 1 1 992h . The SNR of the reduced 
spectra depends on wavelength, and is about 60 and 30 in Ha (6563 
A) and H(5 (4861 A) regions, respectively. 


2.3 Spectral Analysis 

As noted in S ection [T] the S233 ioni zing star was classified 
as B1.5 II by (iHunter & MassevI [l99Qh . We classified this star 
as B0.5 V by comparing our smoothed and degraded to R = 
4000 spectrum with the corr esponding part of the digital atlas of 
IWalborn & FitzpatrickI (ll99Qh (see Figure [T]). Degradation was per¬ 
formed using the GNU Data Language routine ‘REBIN.PRO’. The 
main criteria to estimate the spectral type of the ionizing star were 
approximate equality of the strengths of the Si III 4552 and Si 
IV 4089 lines and the presen ce of the very weak He II 4686 line 
dWalbom & FitzpatrickI [T^ . The ratio of Si III 4552 and He I 
4837 lines was the primary criterion to estimate the luminosity 
class. It seems that the modest disagreement in the spectral clas¬ 
sification is likely due to differences in the resolution and perhaps 
SNR, together with some degree of subjective judgment in the com¬ 
parison by eye with standard spectra. 

To confirm the spectral classification we also performed the 
spectroscopic analysis using stellar atmosphere models from the 
BSTAR2006 grid with the sol ar metallicity and micr oturbulent ve¬ 
locity of Uturb = 2 km s“^ (iLanz & HubenvlE00% . The steps of 
the effective temperature Tgff and surface gravity (log g) in the grid 
are 1000 K and 0.25 dex, respectively. 

At the first step of the spectral analysis the projected rotational 
velocity v sin i was obtained by comparing observed line profiles of 
helium and metal lines with synthetic ones. Synthetic profiles were 
computed with the values of Tgff = 29000 K and log g = 4.0 dex 
which correspond to a B0.5 V star (ISim6n-Diazll2010l) . The com¬ 
parison was performed by eye using synthetic spectra calculated 


by th e ROTIN3 program (provided with SYNSPEC (iHubenv & LanzI 
Il992h ) for three vsini values: 10, 20 and 30 km s The rela¬ 
tively large step of 10 km s“^ allows to easily determine how well 
the observed spectrum is fitted by a synthetic spectrum. Spectra 
with vsini < 10 km s“^ and vsini > 30 km s“^ gave a very 
poor fit. The resulting value of v sin i and its standard deviation are 
13 km s“^ and 5 km s“^, respectively. 

At the second step of the spectral analysis the observed spec¬ 
trum was compared with a number of synthetic spectra computed 
using BSTAR2006 grid for different values of Tgff and iogg. We 
compared observed profiles of Si III AA4552, 4567, 4574, 4716, 
4813, 4819, 4829, 5740 and Si IV AA4089, 4116, 4631 and cores 
and wings of Balmer lines (Ha, H/3, H(5, H7). To characterize how 
well the observed profiles are fitted by model spectra we used the 
quantity 


''^lines f nf \ ^ 

1 ^ / V] ~ Vi obs I 

^lines y J 


( 1 ) 


where rijy is the number of wavelength points in the spectral line i, 
mines is the number of spectral lines used in the analysis, y^ohs is 
the observed flux in the j-ih point of the spectral line z, y] is the 
synthetic flux, Wi is a weight which corresponds to the spectral line 
z, and cr^ = (SNR)“^ accounts for the signal-to-noise ratio of the 
spectral line z. For all lines except Si IV A4089 a value of weights 
equals to 1.0 were used. For the Si IV A4089 line we used a weight 
of 0.5 because it was blended with the Oil line. 

The minimization of provides the maximum likelihood 
estimate of the model parameters, Tgff and iogg. In our case 
the model with Tgff = 28000 K and Iogg = 4.0 exhibit the 
minimum of These model parameters were used to evalu¬ 
ate silicon and helium abun dances based o n the curve-o f -grow th 
method (see, for example, iKilianl (Il992h : ISimon-PfaTl (1201 dri . 
In this method the abundance and Uturb for a given Tgff and 
log g are varied while the same abundance is not obtained for all 
lines. We used a grid of synthetic spectra for various abundances 
and Uturb values that w ere computed with the code SYNSPEC48 
(iHubenv & LanzI IT992h . The number densities of Si in the grid 
are log(Si/H) = —4.75, —4.57, —4.45, —4.3 2, —4.15 dex where 
log(S i/H) = —4.45 is the solar abundance (iGrevesse & Sauvall 
Il998h . The He/H values are varied from 0.05 to 0.15 with 0.025 
step. Using a linear interpolation on the grid with Uturb = 6 km s“^ 
we have obtained the log(Si/H)-EW diagram that is shown in Fig¬ 
ure |2] The log(Si/H)-EW dependence was approximated by the 
straight line that is shown as solid line in Figure |2] The linear ap¬ 
proximation of log(Si/H)-EW is given by log(Si/H) = —4.60 zb 
0.05 + (0.2 zb 0.5)EW, where after a ‘zb’ sign standard errors are 
given. It should be noted that Si III AA4813, 4819, 4829 lines were 
excluded from the Si abundance analysis because this triplet tends 
to give different results. This discrepan cy could be related to th e 
problems of a Si model atom (see e.g. iBecker & Butlej (Il990h ). 
Using the similar diagram, we also obtained the helium abundance 
He/H = 0.10 zb 0.01. 

At the final step of spectral analysis we obtained a x^ dis¬ 
tribution for Uturb = 6 km s“^ and log(Si/H) = —4.60. This 
distribution is shown in Figure [3] According to this distribution, 
we selected all models with a x^ below the critical value. The 
critical value is a minimal x^ value for models on the edge of 
the grid. The resulting values of Tgff = 28000 zb 1100 K and 
log g = 4.15zb0.27 dex are obtained by averaging selected models, 
weighted by exp (—0.5x^), and their standard deviations. 

A comparison of our derived parameters with the calibration 
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Figure 2. Silicon abundance vs. EW diagnostic diagram obtained for 
Uturb = 6 km s“^. The solid line represents the linear regression fit. Un¬ 
certainties in the individual line abundances (propagated from the standard 
deviations in the measured EWs) are indicated as error bars. Because Si IV 
A4089 line was blended with the Oil line, second to last point has large 
uncertainty. 



Te , kK 

Figure 3. difference between the observed spectrum and spectra com¬ 
puted with Uturb = 6 km s“^ and log(Si/H) = —4.60. Abundances of 
other chemical elements are solar. The values are quantified by shades 
of grey. The cross marks the position of the model with minimum value of 


of lCrowth^ (Il997h enables us to suggest that the S233 ionizing star 
is a B0.5 V star. The He abundance also indicates that most likely 
this star is still on the main sequence. We estimated the mass of the 
S233 ionizing star to be M/Mto, = 13 zb 1 by interpolating the stel¬ 
lar e volutionary tr acks of ci are l2004h in the log g-Teff plane (see, 
e.g. jKilia^ (Il992h ) and using the Monte-Carlo error propagation. 
An estimate of the age is much more uncertain with the standard 
deviation equal to the value itself: (3 zb 3)10® yr. Nonetheless, one 
may say that the star is rather young and spent about 0.2 zb 0.2 of 
its lifetime on the main sequence. 

To estimate the distance to the ionizing star we used the K- 
band magnitude K = 9.63 obtained from 2MASS, absolute mag¬ 
nitude Mk = —2.512 teertelli et al.l2009b and interstellar absorp- 
tion Ak — Av^/8.8 — 0. 3 5 acc ording to the reddening law of 
Cardelli. Clavton. & MathisI (Il989l) and Ay = 3.1 estimated by 
Dobashi et al. (12005 ). The estimated distance 2.3 zb 0.4 kpc is in 
agreeme nt with published dist ances to the S231-235 complex. For 
example, IChan & fT^ (Il995l) estimated the distance to the S233 
region as 2.3 zb 0.7 kpc. 

The range of the effective temperature of the ionizing star 
(Teff = 28000 zb 1000 K) limits the range of the number of 
ionizing photons logA^c, to between 47.0 to 47.34. This range 
is obtained by l i near e x trapolation of values on Teff in Table 4 
of iMartins et akl (l2005h . iHunter & MassevI (Il990l) give two esti¬ 
mates of the number of ionizing photons: logiVc = 47.45 and 
log A'c = 47.59, the values are estimated from radio continuum 
and Ka observations, respectively. However, we note that they 
adopted a distance to S233 of 4.7 kpc. If we adopt the distance of 
2.3 kpc, then the values of the number of ionizing photons should 
be scaled by (4.7/2.3)^ = 4.176. This gives log A^c = 46.82 and 
log Nc = 46.97. The value log Nc = 46.97 is almost equal to the 
lower estimate of the number of ionizing photons with the stellar 
parameters obtained by us. This indicates that the ionizing star has 
a sufficient Lyman continuum flux to account for the observed Ha 
intensity and the radio flux density of the S233 H II region. Most 
likely estimation of the number of ionizing photons for the S233 
H II region is log iVc = 47.0. 


2.4 Radial velocity measurements 

The radial velocities were measured for 150 absorption lines us¬ 
ing DECh20. The rad ial velocities were obtained interactively with 
a mirroring method dParimucha & Skod^l2007h . i.e by means of 
the best match of the line profile with i ts mirror. Rest wavele ngths 
for each line were taken from VALD dPiskunov et al.l[l9^ . The 
derived radial velocity of the S233 ionizing star is —17.5 zb 1.4 
km s“^. This velocity is consistent with the velocity of the sur¬ 
rounding molecular cloud (from —15 to —20 km s“^), thus the 
physical association with this molecular cloud may be assumed. 

There are also the nebular lines in the Echelle-spectrum of the 
ionizing star. Using [SII] AA6716, 6731 and [Nil] AA6548, 6583 
lines we estimated the radial velocity —10.7 zb 1.7 km s“^. Be¬ 
cause the nebular lines are used, the derived velocity represents the 
ionized gas in the S233 region. 

Analyzing the velocities of the star, ionized and molecular gas, 
we conclude that the ionizing star is moving together with the par¬ 
ent molecular cloud, but the ionization front is moving slower in 
the direction to the observer with a velocity of 6.8 zb 2.2 km s“^ 
with respect to the ionizing star. This value indicates the expansion 
velocity of the S233 H II region. 

3 RESULTS 

3.1 Nebulosity in the S233 region 

3.1.1 Optical images 

The S233 region is seen as an optical H II region. In the DSS Red 
image the optical nebula exhibits almost a spherical morphology 
(see Figure |4]) which is typically found in H II regions. The DSS 
Red image traces primarily ionized gas. The ionized region does 
not have a clear edge at the north-east while there is a clear edge 
seen in the south-west side. 

The DSS Red image also shows condensation at the west of 
the central ionizing star. This condensation is also seen in near-IR 
and mid-IR images (see Section [3. 1.21 and [3. 1.31 ). The condensa¬ 
tion is also detectable in the DSS Blue image, which cannot be 
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explained as emission directly of the interstellar dust. If we assume 
the emission is the reflected light of the central ionizing star, the 
condensation detected in the DSS Blue image supports that the re¬ 
gion is physically associated with the ionizing star. 


5.7.2 WISE mid-IR images 

Figure [5] shows the mid-IR images of the S233 regio n. These im¬ 
ages were obtained from the all-sky data release dWright et al.l 
bold) of the Wide-held Infrared Survey Explorer (WISE) at 4.6, 
12 and 24 jam. The images reveal the existence of diffuse nebulos¬ 
ity to the south-west of the ionizing star. The main contribution to 
the emission at 12 p.m comes from polycyclic aromatic hydrocar¬ 
bons (PAHs) in a photodissociation region (PDR). Usually a PDR 
appears as a relatively thin region of mid-IR emission along the pe¬ 
riphery of an H II region, where PAH is excited enough to glow, 
but not to be destroyed by UV-radiation 

In our case we do not see a sharp boundary of the PAH emis¬ 
sion of the S233 H II region. This could be due to projection effects 
together with the non-uniform distribution of the ambient material. 
If we overlay the DSS Red image with the WISE 12 p,m image (see 
Eigures|5]and|4]), we hnd that the optical emission of the ionized gas 
overlaps with the PAH emission recognized in the 12 p.m image. In 
the case of uniform distribution of material around a spherically ex¬ 
panding H II region, the bright PAH emission should surround the 
ionized gas. In the present case the PAH emission overlaps with the 
optical emission of the ionized gas, which could be due to projec¬ 
tion effects. This comparison suggests that the H II region interacts 
with the ambient material on the far side from the observer, so the 
region of the PAH emission is relatively extended. Therefore, the 
conclusion is that the ambient material is denser in the south-west 
region: this situation is consistent with the cometary-like structure 
of the ionized hydrogen region, seen in red optical images. 

The nebulosity contains a number of condensations, which are 
labeled on the Pigure|5]as SI, S2 and S3. The bright point source 
below S 1 is likely a foreground star as it coincides with a star seen 
in the optical image, identified as USNO-A2 1200-03587620 and 
2MASS source 05382867-1-3550388 with a magnitude 12.5 in the 
K band. 

In Pigure|5] we see the emission of the ionizing star at 24 p,m. 
In order to check whether this emission is spatially separated from 
the nearest infrared peak SI, we obtained i mages with a higher spa - 
tial resolution using the HiRes technique (iMasci & Eowlejl2009l) . 
The resulting HiRes image of the S233 region shows a clear sepa¬ 
ration between the SI condensation and the ionizing star emission. 
The SI condensation is also seen in the DSS Blue band image as it 
refiects the light of the ionizing star. Prom this fact we can conclude 
that SI is physically associated with the ionizing star in the S233 
region (see Section U.1.11 ) 


3.1.3 UKIDSS near-IR images 

The UKIDSS survey (iLawrence et al.l 1200% is the next genera¬ 
tion near-infrared sky survey, the successor to 2MASS. The sur¬ 
vey instrument is WPCAM (ICasali et al.l l2007h mounted on the 
UK Infrared Telescope (UKIRT) in Hawaii. The high resolution 
of WPCAM (0.4") gives a high quality images of the particular re¬ 
gions of the sky. The S233 region lies in the area covered by the 
Galactic Plane Survey (GPS) of UKIDSS. 

Pigure|4] shows the image in the WPCAM K band along with 
the optical images. It is clearly seen that in the S233 region there 


are well-pronounced structures having forms of bright rims elon¬ 
gated perpendicular to the ionizing star. Within these rim structures 
the condensations SI, S2 and S3 are pronounced. The fine struc¬ 
ture of these condensations is clearly seen in this image. The SI 
condensation has a brightness gradient, revealing the brightest part 
on the ionizing star side and the darker part in the opposite side. 
A number of molecular clouds as sociated with similar structures 
and co ndensations were studied bv [Thompson. Urquhart. & White! 
(l2004 . 

3.1.4 IRAS far-IR data 

Pinally, the IRAS source 0535I-I-3549 is located in the vicinity of 
the S233 region. The angular separation between the IRAS source 
0535I-I-3549 and the SI condensation, discussed in Section [3. 1.21 
is 17". Taking into account the uncertainty of the IRAS source po¬ 
sition (major radius of the error ellipse is 31"), we assume that 
the IRAS source 05351-1-3549 is associated with the SI peak seen 
in the mid-IR images of WISE and near-IR images of UKIDSS. 
This IRAS source has the following fiux densities: F 12 — 1.09 Jy, 
F 25 = 11.5 Jy, Fqq — 184 Jy, F'loo = 414 Jy. It is no¬ 
table that the fiux density increases a s the wavelength increa ses. 
The colours (according to definition of lWalker & Cohen! (!l988!) ) of 
IRAS 05351-^3549 are [12] - [25] = 4.12, [25] - [60] = 4.89. 
These colours correspond to th ose of a Young Stellar Objec t in the 
colour-colour diagram given bv !Straizvs & Laugalvs! (!2007!) . 

One may notice that IRAS 05351-f3549 is marked an extended 
source in the IRAS catalog. This extension is probably due to the 
association of another infrared source with the SI clump. More¬ 
over, the central ionizing star, which is located close to the S1 peak 
{d ~ 24"), emits radiation at mid-IR wavelength and falls into the 
ellipse of IRAS source position error. The angular resolution of the 
IRAS image is not sufficient to resolve the S1 peak and the ionizing 
star and to construct their reasonable spectral energy distributions. 
But mid-IR images from WISE show that these peaks can be dis¬ 
tinguished (see Pigure|5])- 

The Improved Reprocessing of the IRAS Survey (IRIS) im¬ 
ages at 12, 25, 60 and 100 qm show significant brightness and 
extension of the emission of the S233 region, which are compa¬ 
rable to the other regions of star formation within the entire S231- 
235 complex. Ear-IR emission of the S233 region has roughly the 
same rounded morphology with a radius of ~ 4' at all four IRAS 
bands, with the central peak coinciding with the position of IRAS 
05351-^3549. 

3.2 Physical parameters of dust 

In order to estimate the physical parameter s of th e dust we used 
the method described in ICichowolski et^ ] ( I 2001 I 1 . According to 
this method, the total integrated infrared luminosity of IRAS 
05351-1-3549 is as follows (assuming the distance to the S233 re¬ 
gion as 2.3 kpc): 

Lir = 1.58Ftotdgpc = 2.280 x lO" L® (2) 

Where TTot = 1.3(i^i2 + 7 ^ 25 ) + 0.7(7^25 + T^eo) + 0.2(7^6o + 
F'loo) = 272.81 Jy is the integrated source flux density; 
7 ^ 12 , F 25 , Feo, F'loo are the flux densities of the IRAS source at 12, 
25, 60 and 100 qm. 

The calculation of the total mass of radiatively heated dust 
leads to the following results: 

Md = mnF60dkpc(Bn ® - 1) = 0.7 M© (3) 
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Figure 4. UKIDSS near-IR (K-band) and DSS optical (Red and Blue) images of the S233 region. The center of the image 
coincides with the ionizing star at RA(2000) = 5^38™31!5, Dec.(2000) = 35°5lD9''. The image size is 5' x 5'. The 
condensations SI, S2 and S3 discussed in the text (see Section [+11 are marked with labels. 





Figure 5. WISE mid-IR images of the S233 region. The center of the left image coincides with the S233 ionizing star at 
position RA(2000) = 5^38™31!5, Dec. (2000) = 35°51'19". The image size is 5' x 5'. The left image shows distribution 
of 3.4 p.m emission, the two small images at the right shows the distribution of 12 turn a nd 24 p.m emission. The resolution 
of 12 and 24 p.m images was improved with the HiRes technique (iMasci & Fowleill2009h . The condensations SI, S2 and S3 
discussed in the text (see Section l+H are marked with labels. 
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Figure 6. Contours of the (2-1) emission of the S233 region with ve¬ 
locity — 15.6 km s“^, overlaid on UKIDSS 2.2 |am K-band image, showing 
interaction of the S233 H II region with the ambient molecular cloud. Con¬ 
tour levels start from 1.5 K, stepping by 0.5 K. Infrared condensations SI, 
S2 and S3 are marked with labels. 



Figure 7. Line prohles of CO in the direction of S233 (q;j2000 = 
05h38"^28^5 (5 j2ooo = +35°50'51'0. The black solid, gray solid and 
black dashed lines represent the ^^CO (2-1), ^^CO (1-0), and ^^CO (2- 
1) lines. The data cubes are convolved to the same beam size of 47” for 
comparison. Hanning smoothing by 3 points was applied to the ^^CO (1- 
0) line. The prohle of ^^CO (2-1) is multiplied by a factor of 10 to reveal 
low-brightness emission of the clump. 


In this formula n is the emissivity index, related to the absorp¬ 
tion efficiency of the dust (/cv oc normalized to dust ab¬ 
sorption coefficient 40 cm^g“^ at 100 p.m). We used a value of 
n = 1.5 which is typical for the H II regions (iRodon et al.ll2010h 
and mi .5 = 0.3x10“®. Fqq is the fiux density of the IRAS sources 
at 60 p.m and Bn is the modified Planck function, given by the 
equation Bn = 1.667^+^ (Fioo/i^eo) = 22.433. _ 

Using a dust-to-gas ratio of 0.01 dPraine et al.ll2007h we esti¬ 
mated the gas mass to be 70 M©. However, this mass is not neces¬ 
sarily equal to the gas mass that in principle can be estimated from 
CO emission, since infrared emission may trace a different volume 
than that traced by CO emission. 

The last step is a calculation of the dust temperature using the 
following formula: 


Td 


95.94 
In Bn 


30.8 K 


(4) 


This value is typical in star-forming region s dCappa et al.l l2008l : 
ICichowolski et al.ll2009l : IVasauez et al .11201 Oh and consistent with 
the gas temperature in the nearb y star-forming regio n S235 
(iKirsanova et ^1201 4l) and S233-IR (ISchrever et alJI 1996b . 


the S233 H II region (see Figure[6]). It is possible that the ionization 
front moves perpendicular to the observer in this velocity range. It 
should be noted that there is a large amount of material on the west 
of the S233 and there is a void at the east. 

It is possible that the SI peak and IRAS source are the result 
of the interaction between the H II region and ambient molecular 
gas at a velocity of ~ —16 km s“ ^. There also exists a compact red- 
shifted component of gas at a velocity of ~ —13 km s“^, described 
in Section [T4l The SI peak also could be physically associated with 
this component. It may not be possible to draw any firm conclusion 
about whether the S1IR peak is physically related to the molecular 
gas at a velocity of —16 km s“^ or to a redshifted component of 
gas at a velocity of —13 km s“^. 

The line intensity ratio of ^2CO(2-l)/^^CO(2-l) in the S233 
region ranges from 3 to 15, depending on the velocity shift 
from the line peak. Using the main beam antenna temperature of 
the line peak, we obtained a lower limit for the line ra¬ 

tio ^^CO/^^CO ~ 3. Using the values of antenna temperature 
away from the line peak, we got the higher line intensity ratios 
^^CO/^^CO ~ 10 — 15. This is direct evidence that the line 
is optically thick. 


3.3 Ambient molecular cloud 

The S231-235 complex is a warm giant molecular cloud with 
brightness temperatures of the (2-1) line up to 50 K. Most 
of CO emission has main beam antenna temperatures Tmb ^ 16 K. 
The compl ex was observed ea rlier in the (1-0) and (1- 
0) lines bv iHever et al.l (Il996h . 

_ ^It appears that molecular gas in the S233 region (iHever et al.l 

Il996h is the part of the large-scale filament that lies between 
the expanding S231 H II region and the void at the west of 
S233. This void shows a circle - like m orphology in maps given by 
IPame. Hartmann. & Thaddeui (l200lh . 

A hint suggesting a physical interaction between the S233 H II 
region and ambient molecular clouds is seen in the velocity range 
from —16.1 to —14.5 km s“ ^. The irregular border of CO emission 
at the west of the S233 region may be physically associated with 


3.4 Molecular clump in the S233 region 

Figure |7] shows the profiles of the CO molecular line in the di¬ 
rection of the S233 region (aj 2 ooo = 05^38'^28®5; (5j2ooo = 
+35°50^51^0- The primary component exhibits velocities in the 
range from —20 to —15 km s“^, and a weaker red-shifted compo¬ 
nent is peaked at a velocity of ~ —13 km s“^. The red-shifted com¬ 
ponent is pronounced in the ^^CO (2-1), ^^CO (2-1), and ^^CO (1- 
0) lines. Channel maps of the red-shifted component are presented 
in Figure [8] 

From the integrated image of the red-shifted component we 
estimate the spatial extent of this component using Gaussian fits to 
the spatial intensity profile along two perpendicular directions (see 
Figure |9]). The source size is defined as FWHM of the fitted Gaus¬ 
sian functions. Derived values are di = I'll" ± 4" and d 2 = 
1'4" ± 6.6", where the first intensity profile is obtained along 
the line passing through the ionizing star and the intensity peak 
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Figure 8. (2-1) channel maps of the clump in the S233 region. The velocity of each panel is shown at the top-left 

corner. The ellipse represents the visual extent of the S233 H II region. The cross indicates the position of the ionizing star. 
The vertical and horizontal lines exhibit the position of the intensity peak of the clump at —13.18 km s“^. 


of the red-shifted component. The second profile d 2 is obtained in 
the perpendicular direction. The source size appears to be approxi¬ 
mately twice that of the telescope beam size 0 = 32", therefore the 
red-shifted component is not a point source. We conclude that this 
red-shifted component corresponds to a separate clump of gas. As¬ 
suming the distance to the S233 region D = 2.3 zb 0.3 kpc, the di¬ 
ameter of the clump is tii = 0.26zb0.037pc, d 2 = 0.23zb0.038 pc. 

The main beam brightness temperatures of the clump emission 
are 3.42 K, 2.79 K, and 0.26 K in the (2-1), (1-0), and 

(2-1) lines, respectively. The clump emission shows a clear 
velocity gradient (see Figure [8]). The minimum velocity at which 
the clump becomes visible in the (2-1) line is —14.3 km s“^, 
and the maximum one is —11.5 km s“^. The position of the in¬ 
tensity peak of the clump shifts by ~ 1' when the radial veloc¬ 
ity changes from from —14.3 km s“^ to —11.5 km s“^. At the 
distance of 2.3 kpc this velocity shift corresponds to a gradient of 
u ~ 1.8 km pc“^. 

The clump emission is located around the SI IR condensa¬ 
tion described in Section 13.1.21 The existence of the gradient in 
the clump emission is probably a result of the interaction with the 
expanding H II region. The position of the clump at the median ve¬ 
locity (—13 km s“^) coincides with the position of the SI mid-IR 
peak. As the velocity increases to — 11 km s“^ the intensity peak of 
the clump shifts to the west relative to the S1 condensation and the 
edge of the H II region. As the velocity decreases to —14.6 km s“^ 
the intensity peak of the clump shifts to the south relative to SI 
(see Figure [8]). These trends show that the clump has a complex 
morphology with the S1 mid-IR peak in the center. 

The clump is seen in three CO lines, but the emission is rather 



Figure 9. Contours of the (2-1) emission, integrated from —14.6 

to —11.0 km s“^, showing the molecular clump in the S233 region. The 
contour levels span from 2.22 K km s“^ to 6.66 K km s“^. The background 
is the WISE 12 pm image. Filled ellipse shows uncertainty in position of 
IRAS source 05351+3549. Intensity profiles along the perpendicular dotted 
lines were used for the clump size determination (see Section lT^ . 

weak, so the noise levels prevent a detailed analysis of CO excita¬ 
tion using three lines simultaneously. So we applied only a simple 
analysis assuming LTE conditions and a uniform single gas layer. 
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Table 1. Parameters of CO line profiles of the S233 molecular clump in 
direction aj 2 ooo = 05*'38"’28?5, 5j2000 = +35°50'51". The first col- 
umn is the observed main beam temperature of line, the second column is 
the central velocity of observed line and the third column is the FWHM 
of observed line. Standard error for each parameter is obtained from the 
Gaussian fit. 


CO line 

^peak 

(K) 

Vo 

(kms-i) 

Al/ 

(kms-i) 

i^CO (1-0) 

2.79 (0.6) 

-13.48 (0.154) 

2.16(0.39) 

12C0 (2-1) 

3.42 (0.07) 

-13.21 (0.015) 

2.20 (0.04) 

(2-1) 

0.26 (0.04) 

-13.36 (0.081) 

1.88 (0.20) 


These assumptions allow us to estimate an optical depth using the 
main-beam brightness temperatures of the and in (2-1) 
lines (see Table [T]) using the following equation: 


Ti 2 _ 1 — e 

7 ^ ~ 1 


(5) 



Observer 


Figure 10. Schematic representation of the S233 H II region constructed 
from the UKIDSS image (Figure IdJ and velocities of the ionizing star, ion¬ 
ized gas (Section ES and molecular clump (Section EIJ. Dotted region 
shows molecular gas with the range of velocities from —20 to —15 km s“ ^. 


here X = 60 zb 10 is th e relative abundance of [^^CO]/[^^CO] 
at the distance of 2.3 kpc (iLanger & Penziaslll990h . Using the ob¬ 
served value of T 12 /T 13 = 13.15 zb 2.04 and solving Equation ([5]) 
numerically, we estimate that the optical depth of the (2-1) 

line is in the range from 3.3 to 6.6. This shows that the line 
is optically thick. 

For the extended sources the sum of the brightness tempera¬ 
tures of the optically thick line and background, T 12 + Tbg, pro¬ 
vides an estimate of the gas kinetic temperature. In the case of 
the molecular clump both the (2-1) and (1-0) lines have 

rather low values of main beam temperature Tmb — 3 K. Cor¬ 
recting the observed radiation temperature for the Planck law, i.e., 
Tmh = — (hiy/k)/[e:>cp(hiy/kT) — 1], the derived excita¬ 

tion temperature for the (2-1) line is 8 K for an observed peak 
line intensity of 3.4 K and assuming only the cosmic background 
of 2.73 K. A value of 8 K is a lower limit to the gas kinetic tem¬ 
perature, because (1) the transition may be sub-thermally excited; 
(2) the continuum background may include some thermal dust con¬ 
tinuum; and (3) the source may be only partially resolved so the 
observed brightness temperature is beam-diluted. The latter is in 
accord with recent observational results showing that the formation 
of low mass stars tend s to occur in filaments w ith a characteris¬ 
tic size of about 0.1 pc (I Arzoumanian et al]l201lh which is smaller 
than the linear size of the beam in the present observations. Obser¬ 
vations with a higher spatial resolution and signal-to-noise ratio are 
necessary to rev eal the situation 

Following iMacLaren et all (Il988h . we estimated the virial 
mass of the clump using the linewidth of the line (1.88 

km s“^) and the mean clump radius of 0.12 pc. Derived values 
range from 53 to 89 M©, depending on the assumed density distri¬ 
bution as a function of the radial distance from the clump center r. 
The virial mass is calculated to be 53, 80 and 89 M© in the cases 
of a constant density, density oc 1/r and oc 1/r^, respectively. All 
these values are of the same magnitude as the mass of the IRAS 
source 05351-1-3549, derived from IRAS flux densities (70 M©, see 
Section \3^ for details). This agreement indicates that the clump 
may be gravitationally bound if the IRAS source is physically as¬ 
sociated with the CO clump. 


4 DISCUSSION 

4.1 Morphology of the S233 H II region 

It is seen in the optical images (see Figure 0]) that the S233 region 
has a co metary shape w ith the sharp edge situated in the south-west 
(see, e.g. lCohenl(ll98Cll) ). This morphology could be due to the mo¬ 
tion of the ionizing star relative to the ambient gas, or to inhomo¬ 
geneity of the ambient cloud. In the Section [T4l we have shown that 
the velocity of the ionizing star (Ulsr = —17.5 zb 1.4 km s“^) lies 
within the velocity range of the ambient molecular gas (Ulsr from 
—20 to —15 km s“^) that surrounds the star. So, it is more plausible 
that the central ionizing star was born in the molecular cloud, not 
that the star entered this cloud from the outside. Thus, the cometary 
morphology the S233 region seen in DSS Red image probably re¬ 
sults from an inhomogeneous distribution of the ambient molecular 
gas (see Figure[T0l). If the material to the south-west from the star is 
denser than the material to the north-east, then the photoionization 
and shock waves propagating to the south-west meet more resis¬ 
tance, thus producing the brighter part of the cometary-Iike nebula. 

The condensations SI, S2 and S3 discussed in the Section 
13. 1.21 might be formed by the expansion of the H II region. In the 
vicinity of S1 the gas clump could be dense enough that its material 
cannot be ionized or pushed away by expansion of the S233 H II re¬ 
gion. The images clearly show that the condensations are oriented 
tangentially with respect to the direction towards the ionizing star. 
Thus, these condensations might be created by the photoionization 
and shock waves expanding from the S233 H II region. 

4.2 A scenario for star formation in the S233 region 

In the S233 region we see evidence for interaction between the H II 
and the ambient molecular cloud, which could trigger star forma¬ 
tion. There is a compact infrared cometary-Iike condensation SI, 
which faces the ionizing star and exhibits emission at 24 qm which 
comes primarily from heated dust. This condensation coincides in 
position with the IRAS source 05351-1-3549. This source displays 
flux density increasing with wavelength, with colours correspond- 
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Table 2. Characteristics of the one-dimensional model simulating the 
‘collect-and-collapse’ process for the S233 H II region. Parameters used 
for calculation of the triggering process are: logA^c = 47.0, -Rgheii = 
0.67 pc, Tkin = 25 K (see explanations in the text). 


riQ, cm ^ 

102 

5 • 102 

103 

3 • 103 

4- lO-i 

i^s(0),pc 

0.74 

0.25 

0.15 

0.05 

0.01 

-^frag’ 

15.9 

6.6 

4.5 

1.9 

0.6 

^frag? Myr 

9.1 

4.4 

3.2 

1.5 

0.6 

Mfrag, M© 

560 

270 

200 

100 

40 

Age, Myr 

0.03 

0.06 

0.09 

0.2 

0.6 


ing to those of a YSO (see Section [5X4]). There is also a molecular 
clump in the direction of the SI condensation, seen in emis¬ 
sion red-shifted relative to that of the parent molecular cloud. 

It is known that the expansion of an H II region can trigger 
star formation. There are different scenarios of star formation trig¬ 
gered by the H II region expansion. One of t he commonly con¬ 
sidere d scenarios is the ‘collect-and-collapse’ (lElmegreen & Lada! 
I 1977 I) . In this scenario neutral gas in the layer between the ioniza¬ 
tion and shock fronts becomes gravitationally unstable which even¬ 
tually leads to star formation. In the ideally symmetric case this 
neutral gas is expected to be observed as a ring structure around the 
H II region. The young stellar clusters are expected to be located 
between the P DR and the parent molec ular gas. 

Following [Kirsanova et al.l ([20141) we consider the S233 H II 
region within the framework of the theor etical one-dimensional 
model described in [Whitworth et al.[ ([l994l) . The model allows an 
estimate of the main characteristics of the ‘collect-and-collapse’ 
process: starting time of the dense shell fragmentation (tfrag), the 
radius of the shell (i^frag) at this moment and the mass of the frag¬ 
ments (Mfrag). We also calculate the radi us of the initial Stromgren 
sphere i?s(0) ([Pyson & Williamsl[l997[) and the H II region age 
(using Eq.4 from I Whitworth et ^ ([l994[) ). However, the S233 re¬ 
gion exhibits a non-uniform distribution of the ambient molecular 
material. Thus, the one-dimensional model with a homogeneous 
ambient density distribution cannot be considered as strictly appli¬ 
cable to this region. Since there is no dense gas to the north-east 
of the H II, the pressure in the H II region is lower than in the 
model. Therefore we use the model of [Whitworth et al.[ ([ 19941) to 
make only rough estimates. 

Input parameters of the one-dimensional model are the num¬ 
ber of ionizing photons, the size of the H II region, the kinetic 
temperature, and the gas density. The number of ionizing photons 
for the S233 ionizing star is Nc — 10^^ ° (see Section [23]) . For 
the calculations we used a fixed value of the H II region radius. 
The angular diameter of the S233 reg ion is 120", taken from the 
optical measurements ([SharDless[[l959[) . Publishe d distances to the 
objec ts in this complex vary from 1.6 to 2.3 kpc ([Reipurth & Yan[ 
[2008[) . Thus, the linear radius of the S233 region vary from 0.46 pc 
(D = 1.6 kpc) to 0.67 pc {D — 2.3 kpc). We adopt the value of 
0.67 pc as a representative radius of the S233 H II region. We adopt 
25 K as a typical te mperature in the S231-235 complex accordin g 
to measurements of [Kirsanova et al.[ ([2014[) : [^chrever et~^ ([l996[) . 
Note that the parameter estimates do not depend critically on the 
adopted temperature. For the gas density we consider several val¬ 
ues from 100 cm“^ to 40000 cm“^. The age of the H II region 
becomes comparable to tfrag only with the highest assumed value 
of the gas density. Results of the calculations are given in Table [2] 

Analyzing the results of these calculations, we conclude that 
the ‘collect-and-collapse’ scenario in the S233 region could work 


if the average gas density in the region has the value no ^ 4 • 
10^ cm“^. With the density this high the fragmentation time be¬ 
comes comparable to the H II region age. The estimate of the mass 
of the clump in that case is Mfrag ~ 40 M©, which is compara¬ 
ble with the gas mass estimate from the dust emission (70 M©, see 
Section [ 33 ]) . We cannot estimate the gas density in the S233 re¬ 
gion from our data alone. However, the nearby S235 H II region 
exhibits an average ga s density of 500 cm“^ only (see Table 1 in 
[Kirsanova et al.[ (12014[) ). Thus, the value of 4-10^ cm“^ is probably 
too high for the average gas density in the S233 H II region. Based 
on cal culations using one-dimensional model of [Whitworth et al] 
(Il994h . we conclude that the ‘collect-and-collapse’ process is not 
likely to take place in the S233 region. 

Another scenario for triggered star formation is contrac¬ 
tion of a pre-existing clump due to radiation-driven implo- 
sion J Lefiyh & Lazaref j [1994 [Kessel-Devnet & Burkerd [2()03[ : 


[Miao et ^ 


2006 ) or to a shock wave from the expanding H II re¬ 


gion ([Boss[ll995[) . This scenario predicts a random distribution of 


molecular clumps around the H II region. The signature of this 
process is the presence of the bright rims surrounding cometary 
globules and condensations similar t o those considered, e.g. by 
I Thompson. Urquhart, & Whit3 ([20041) . Star formation could take 
place in these globules. In the case of the S233 H II region we ob¬ 
serve a bright rim with a cometary-like infrared condensation SI, 
faci ng the S233 i onizin g star (see Figure[4]). 

iMiao et akl ([2006[) conducted a three-dimensional hydrody- 
namical modeling to study the radiation-driven implosion effects of 
massive stars on surrounding molecular clouds. They showed that 
different morphologies of the bright rim represent different stages 
in the evolution of the cloud. The S1 condensation has a shape that 
is very similar to that of a molecular clum p at stage I, accord ing to 
the classification of evolutionary stages bv [Miao et al.[ ([2006[) . The 
distribution of infrared emission in the WISE wavelength bands 
shows that the gas in S1 is hotter on the side that faces the ionizing 
star (figure [ 5 ]). At present the SI condensation is elongated tangen¬ 
tially but it is possible that eventually condensation SI will become 
more elongated in the radial direction pd w ill become a classical 
cometary globule ([Lefioch & Lazaref3[ 19941) . The morphology of 
the optical and infrared emission and their spectral characteristics 
show that S1 represents a site of active interaction of the ionizing 
star emission with the dense gas clump. The SI peak is associated 
with the IRAS point source, and probably is situated at the isolated 
molecular clump. It is possible, therefore, that star formation in the 
S1 region takes place because the shock wave from the H II region 
has compressed the clump of ambient material. 


5 CONCLUSIONS 

We have studied a star-forming region which includes an isolated 
B0.5 V ionizing star surrounded by the H II region S233. The 
ionizing star has a sufficiently large Lyman continuum flux to ac¬ 
count for the observed radio flux density of the S233 H II region. 
The velocity of the ionizing star is determined to be Vlsr = 
— 17.5 zb 1.4 km s“^. Comparing the velocity of the ionizing star 
with that of the ambient molecular gas (from —20 to —15 km s“^), 
we conclude that this star was formed from the material of the am¬ 
bient gas cloud. 

In the Section [2!4] we have shown that the S233 H II region ex¬ 
pands with a velocity of 6.8zb2.2 km s“^ and interacts dynamically 
with the ambient molecular cloud material. We find evidence for 
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interaction between S233 and the ambient molecular cloud, which 
could trigger star formation: 

(i) UKIDSS and WISE images reveal an extended region of in¬ 
frared emission containing the compact condensation SI (up to 24 
|am) and the bright-rimmed structure, perpendicular to the ionizing 
star. The bright-rimmed structures with infrared condensations are 
com mon in star formation triggered by r adiatively-driven implo¬ 
sion (iThompson. Urguhart. & Whitel2004h . The dust component of 
the infrared condensation SI reflects visible light (DSS Blue) from 
the ionizing star, which supports the idea of physical proximity of 
the star and the infrared source. 

(ii) IRAS source 05351-1-3549 coincides in position with the 
compact infrared source SI in the WISE and UKIDSS images. The 
flux of the IRAS source increases with wavelength (up to 414 Jy at 
100 jam) and the colours of the IRAS source c orrespond to those of 
a YSO in the colour-colour diagram given bv IStraizvs & Laugalvsl 
(l2007h . The dust temperature of the infrared source, Td = 30.8 K is 
in good agreement with the gas temper ature of star forniing cl umps 
in the nearby star-forming region S235 (iKirsanova et al .l2014h . The 
gas mass of the IRAS source is found to be 70 M©. 

(iii) CO study: The infrared condensation SI coincides in po¬ 
sition with the clump of molecular gas having a mean velocity 
— 13.1 km s“^. The main gas component in the region has veloc¬ 
ity in the interval from —20 km s“^ to —15 km s“^, so the clump 
may have been ‘disturbed’ by the H II region and is now mov¬ 
ing away from us with a velocity which exceeds the velocity of 
‘undisturbed’ gas. The clump emission is observed at velocities 
from —13.4 km to —11.8 km and shows a velocity gradi¬ 
ent (jj ^ 1.8 km s~^ pc“^. The optical depth of the clump in the 

lines is in the range r = 3.3 — 6.6. The virial mass of the 
clump ranges from 53 to 89 M©. 

We investigate the nature of the triggering mechanism for star 
formation in the S233 H II region vicinity. Analysis of calculations 
using the one-dimensional model given in I Whitworth et all (Il994h 
indicates that the ‘collect-and-collapse’ scenario in the S233 region 
could be realized if the average gas density in the region has the 
value no ^ 4-10^ cm“^. This value is probably too high for the 
average gas density in the S233 H II region. Thus, we conclude that 
the ‘collect-and-collapse’ process is not likely to take place in the 
S233 region. 

A more plausible scenario for triggering of star formation in 
the S233 H II region is the ‘collapse of the pre-existing clump’. 
The signature of this process is the presence of the bright rims 
surrounding cometary globules and condensations. In the case of 
the S233 H II region vicinity we observe a bright rim around 
the cometary-like condensation S1 seen in the infrared continuum 
emission. This rim is facing the S233 ionizing star (see EigurelH). 
Therefore, we suggest that star formation of the ‘collapse of the 
pre-existing clump’ type is taking place in this condensation. 
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